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Abstract 

The X-ray properties of a sample of high redshift (z>0.6), massive clusters observed with 
XMM -Newton and Chandra are described, including two exceptional systems. One, at 
z=0.89, has an X-ray temperature of T=ll .5^11 keV (the highest temperature of any cluster 
known at z>0.6), an estimated mass of «(1.4±0.2)xl0'^ Mq , and appears relaxed. The 
other, at z=0.83, has at least three sub-clumps, probably in the process of merging, and may 
also show signs of faint filamentary structure at large radii, observed in X-rays. In general 
there is a mix of X-ray morphologies, from those clusters which appear relaxed and con- 
taining little substructure to some highly non-virialized and probably merging systems. The 
X-ray gas metallicities and gas mass fractions of the relaxed systems are similar to those 
of low redshift clusters of the same temperature, suggesting that the gas was in place, and 
containing its metals, by z«0.8. The evolution of the mass-temperature relation may be con- 
sistent with no evolution or with the "late formation" assumption. The effect of point source 
contamination in the ROSAT survey from which these clusters were selected is estimated, 
and the implications for the ROSAT X-ray luminosity function discussed. 



1.1 Introduction 

Massive clusters of galaxies are rare objects, forming from the high-sigma tail of 
the cosmological density distribution. Their properties are powerful probes of cosmology, 
and give insight into the process of structure formation on large scales. For example, com- 
pared to z=0, and in a ACDM Universe, the number density of virialized halos of M > 5x10'"* 
Mq is predicted to be a factor of fslO lower at z=0.5, or a factor of wlOO lower at z=l, from 
both simulations and the Press-Schechter approximation (eg. Bode et al. 2001). Such a 
large change in number density implies that a significant fraction of the progenitor systems 
at high redshift should be unvirialized, perhaps containing lower mass systems in the process 
of merging. 

The number of genuinely massive clusters known at high redshifts is very small, and 
the number of relaxed high redshift massive clusters is even smaller. Thus the fraction 
which are relaxed, giving direct information on the epoch and mode of assembly of massive 
clusters, is poorly known. Relaxed, massive clusters at high redshifts also offer the best 
opportunities for deriving total masses and gas mass fractions from the X-ray data without 
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too many uncertainties due to hydrostatic equilibrium assumptions, and also with reasonable 
signal-to-noise. 

Detailed analyses of 3 of the clusters presented here can be found in Maughan et al. (2003a,b). 
We use n„,=Q.3, rt\=0.1 and Ho=70 km s"' unless otherwise stated. 



1.2 The cluster sample 

The clusters were discovered in a ROSAT serendipitous X-ray survey (WARPS - 
Scharf et al. 1997, Jones et al. 1998, Perlman et al. 2002). Of 16 clusters at z>0.6 which form 
a complete. X-ray selected sample, 10 have been observed with XMM and/or Chandra. For 
3 of these, the XMM observations were of low quality due to periods of high background. 
The remaining observations, however, represent a major step forward compared to previous 
results from ROSAT (see also Rosati et al. , these proceedings and references therein). 

1.3 Analysis of the high redshift clusters 

Standard analysis techniques were applied to the XMM and Chandra data, taking 
into account the ACIS contamination affecting the low energy efficiency, and using a lo- 
cal background subtraction method to obtain spectra. Images were exposure corrected and 
adaptively smoothed so that all features are significant at the 99% confidence level. 

In general a single temperature was measured (but see CLJ1226.9H-3332 below for a tem- 
perature profile). Two-dimensional /^-profiles were fit to the images, including PSF blurring, 
to obtain values for the core radii and (3. For the clusters which appeared relaxed, and as- 
suming spherical symmetry, the gas mass was derived from an integral of the density profile 
given by the /3-profile. The total gravitating mass was derived by additionally assuming 
isothermality and hydrostatic equilibrium. The total mass within a radius r is then given by 



M«r) = 1.13x10'^/?-^ / ^'j^f Mr,. (1.1) 

keV Mpc 1 + (_r/rc) 

To measure the virial radius, we used eqn (1.1) to find the radius within which the total 
overdensity was 200 times the critical density at the redshift of each cluster. The quoted 
fluxes, luminosities and masses were measured via an extrapolation from the typical detec- 
tion radius of ~0. 3-0.45 r,,„. to the virial radius. 



1.4 Descriptions of the clusters 

7.4.7 A massive, relaxed cluster at z=0.89 

The discovery of CLJ1226.9H-3332 is described by Ebehng et al. (2001). The XMM 
image of CLJ1226.9H-3332 is shown in Fig ll.ll overlaid on a deep I band image. The X-ray 
morphology appears generally relaxed, in contrast to cluster MS 1054-0321, the other high 
temperature cluster known at z>0.8 (Jeltema et al. 2001). The XMM temperature of 1 1 .S^.^'l, 
keV is consistent with the velocity dispersion of Ril 100km s"', based on 15 galaxy redshifts. 
The bolometric X-ray luminosity is 5.4x10'*'' erg s~^ We derive a mass of (1.4±0.2)xl0'^ 
M© within the virial radius. Our analysis of a short Chandra observation confirms the lack 
of significant point-source contamination, and the temperature and luminosity, albeit with 
lower precision (see also Cagnoni et al. 2001). 

A temperature profile (Fig ll.H is consistent with the cluster being isothermal out to 45% 
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Fig. 1.1. Top - XMM contours (0.3-8 keV) overlayed on a Subaru I-band image of 
cluster C1J1226. 9+3332 (z=0.89). The contours are logarithmically spaced. Bot- 
tom - XMM X-ray temperature profile of C1J1226. 9+3332. The solid line is at the 
best fitting global temperature, with la errors shown by the dashed lines. 



of the virial radius. The metal abundance of Z=Q.33'!^^ \qZq and gas mass fraction of 8.6^1,',% 
are consistent with those of local clusters of the same temperature (see below). 

This cluster is unique in being at high redshift, yet massive and generally relaxed. It must 
have been assembled when the age of the Universe was significantly less than 6 Gyr Further 
details of the XMM analysis and the cosmological implications are given in Maughan et 
al. (2003b). 

1.4.2 A massive cluster in formation at z=0.83 

In stark contrast is the highly unrelaxed, massive cluster C1J0152. 7-1357, described 
in Ebehng et al. (2000), Delia Ceca et al. (2000) and Maughan et al. (2003a). The Chandra 
image and a recently acquired deep XMM image are shown in Fig ll.2l Based on the Chandra 
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data, we find that the cluster consists of two major sub-clumps at very similar redshifts 
and of temperatures S.S;!;!} g keV and 5.2^}^ij keV. The total mass of each subcluster is w(5- 
6)xl0^'* Mq . A dynamical analysis of the system shows that the subclusters are likely to be 
gravitationally bound. When merged, the system mass will be similar to that of the Coma 
cluster. 

The deeper XMM image shows lower surface brightness features. To the east is a low 
luminosity system also at the cluster redshift (Demarco et al. 2003, in preparation). To the 
NW is a possible filamentary structure with a point source and one or two low luminosity 
extended sources, probably groups, embedded within it. This structure falls within a ra- 
dius of 1.6 Mpc (the dashed circle), the estimated virial radius of the final merged system, 
suggesting that it is part of the same formation process. The faintest filamentary emission 
is significant at >5a in the soft 0.3-1 keV band and has a hardness ratio consistent with a 
temperature of 1±0.7 keV. We may be observing the warm baryons predicted to lie within 
filaments, based on simulations of structure formation. 

1.4.3 Other clusters at z>0.6 

Chandra and XMM contours of the remaining 5 clusters are shown in Figs ll.3lll.4l 
and ll.5l in order of the most relaxed first. There is a wide range of morphologies, and some 
(eg C1J1559) clearly have bright point X-ray sources nearby. 

An additional cluster, C1J1227.3H-3333 (z=0.766) is shown in Fig ll.6l This is an extremely 
unrelaxed system. Although not in the original ROSAT sample (it fell below a window 
support shadow), it is detected in an XMM image and is a factor ^^12 less luminous than 
C1J1226.9H-3332. There is no clear optical or X-ray core, but rather a system of Ri4 X- 
ray subclumps separated by ^0.5 Mpc which are at a very similar redshift. There is a red 
sequence in the galaxy colours. This is probably a highly non-virialized system seem in 
formation, yet it is detected in X-rays. 

1.5 Metallicities and gas mass fractions 

Most of the metals in a cluster are in the X-ray gas, and so the evolution of the 
metallicity of the gas provides important information on the chemical history of the Uni- 
verse. In Fig ll.Tl we show that our measurements, and those of others, of cluster metallici- 
ties up to z=0.9 are consistent with the canonical value of 0.3 times the solar value found at 
lower redshifts by Mushotzky & Loewenstein (1997). 

The gas mass fractions we measure, within a radius of 0.3r,,,> (where they are the most 
reliable) are consistent with those found by Sanderson et al. (2002) for local clusters of the 
same temperature. Thus the gas, and the metals it contains, were in place by zwO.8 in at 
least some massive clusters. 

1.6 The evolution of the mass-temperature relation 

In FigQ we show the mass-temperature relation for those high redshift clusters 
which appear reasonably relaxed (shown as red squares). The circles are the low redshift data 
points of Sanderson et al. (2002). In both cases the mass is that measured within the virial 
radius, r2oo, and in both studies r2()o is measured from the derived mass profile as the radius 
within which the density is 200 times the critical density at the redshift of observation. The 
low and high redshift samples have different treatments of the temperature profile, however 
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Fig. 1.2. Top - Chandra contours of C1J0152.7-1357 (z=0.833). Bottom - XMM 
greyscale image of C1J0152. 7-1357. The circle denotes the estimated virial radius 
of the final, merged, system. Note the group to the east (at the same redshift as the 
cluster) and the filamentary structure, with embedded sources, to the NW. 



For the high redshift sample we (necessarily) assume isothermality. For the low redshift 
sample, the temperature profile is measured and thus a more accurate mass is derived. 

The upper dashed line shows the best fit to the low redshift points. The lower (red) dashed 
line shows the expected change in normalisation assuming clusters observed at z=0.8 were 
formed at z=0.8, using the relation E{z)M ex T'^l^, where E{z) = (1 +z) ^/(l +z^M+^Ali} +z)'^-^a)- 
It is interesting to note that all the high redshift clusters except the most massive one are con- 
sistent with the low redshift relation, rather than the z=0.8 prediction. It is however the most 
massive cluster (CU1226) which has evidence of isothermality and the most rehable mass. 
The masses of the other high redshift clusters would need to be reduced by w40% in order 
to be consistent with the z=0.8 prediction. Sanderson et al. have shown that an incorrect 
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Fig. 1.3. Left - Chandra contours of ClJl 113.1-2615. Right - XMM contours of 
CU0046.3+8531 
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Fig. 1.4. Left - XMM contours of CU1008.7+5342 Right - XMM contours of 
CU1559.1+6353. 
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Fig. L5. XMM contours of CU 1342.8+4028. 
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Fig. 1.6. XMM contours of C1J1227. 3+3333 (z=0.766), which is adjacent to 
C1J1226.9+3332 on the sky, ahhough at a different redshift. The contours are Un- 
early spaced and the underlying image is in the R band, from Subaru. 



assumption of isothermality when clusters are not isothermal produces, on average, masses 
which are too high by 30% at r2oo- This would explain most of the offset we find. Further 
work is in progress, and will be reported in future papers, but more accurate masses at high 
redshifts, based on temperature profiles, may be required before firm conclusions can be 
reached. 

1.7 Point source contamination of ROSAT luminosities 

The ROSAT detections were obtained at ROSAT off-axis angles of up to 15 arcmin, 
where the PSF degrades to w50 arcsec (FWHM). Thus some point-source contamination 
is expected in all the ROSAT cluster surveys. The degree of contamination depends on 
the signal-to-noise ratio of the detection and the algorithm used for source detection and 
characterisation. 

Based on a preliminary analysis, 2 or 3 out of 9 clusters at z>0.6 in the WARPS survey 
were artificially boosted above the survey flux limit by significant point source contamina- 
tion. Removing these sources from the sample, and using the updated XMM and Chandra 
fluxes, we recomputed the X-ray luminosity function at 0.6 < z <1.1. The resultant shifts 
to lower space densities were within the Poissonian error bars, suggesting that point source 
contamination can be an important effect for individual ROSAT measurements, but the pre- 
vious broad conclusions regarding cluster evolution are unaffected. 
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Fig. 1.7. Cluster metallicity vs. redshift. Circles at low redshifts are from 
Mushotzky & Loewenstein (1997), with 90% confidence error bars, filled (red) 
squares represent this work, with 68% confidence error bars, and the open square 
is from Jeltema et al. (2001). 
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Fig. 1.8. Total mass within r2oo vs. X-ray temperature. Squares (red) represent the 
relaxed clusters in this work at high redshifts. Open circles are from Sanderson et 
al. (2002) at low redshifts. See text for details of the measurement methods and the 
z=0.8 prediction (lower dashed line). 
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